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ABSTRACT 


In the present work, a scheme of geometrically mode- 
ling superquadric solid objects using the constructive solid 
geometry representation has been designed and implemented. A 
superquadric surface enables the designer to model objects 
having fillet- radii and rounded edges using a single bi-para- 
metric equation. Also, a variety of different types of objects 
can be defined by varying the value of the exponent in the 
bi-parametric equatiop. Algorithms for the boolean operations 
of Solids defined by superquadric surfaces have been developed 
using the constructive solid geometry approach. 

In order to render the objects in a realistic manner, 
the ray tracing technique has been used. Also, the shading 
of surfaces has been achieved by using Phong's shading model. 
The entire scheme of geometric modeling has been implemented 
in the form of a menu-driven program on ND— 560 system using 
Tektronix- 41 09 graphics display system. Illustrative displays 
of primitives and objects are presented. 



CHAPTER 1 


INTRODUCTION 


1 . 1 Geometric Modeling - A Concept 

Computer graphics has grown as an area of computer 
applications enormously. It has developed applications in 
every possible field ranging from arts having use in computer 
painting and animation, in medical field for CAT scanning. 
Also, computer graphics has got wide engineering applications 
in the field of Robotics, CAD, CAM, finite element analysis 
and geometric modeling to name a few. 

Geometric modeler is essentially a computer based 
system which provides an efficient tool for creating, editing 
and storing the object representation needed for different 
engineering applications. Display generation is an integral 
part of the geometric modeler. An unambiguous display of 
the solid is an information bank in itself. Majority of the 
processes which use the modelled geometry to compute the 
properties of the object (such as mass, moment of inertia, 
volume etc. ) are application subsystems outside the geometric 
modeler. 

Representations in a geometric modeler are through 
an input language® It may range from simple interactive 
graphics command language to an elaborated block structured 
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language similar to conventional programming languages. The 
sequence of the input language statements is itself an object 
representation . 

1.2 Geometric Definition of Solid 

The core of the geometric modeler consists of the 
symbol structures (representations) designating "Abstract 
solids" which are manipulated to model the desired solid. 

The abstract solid has following characteristics [l]: 

i) RIGIDITY t It should have an invariant configura- 
tion or shape which is independent of solids loca- 
tion or orientation. 

ii) homogeneous THREE DIMENSIONALITY : A solid must 

have an interior. It cannot have isolated or 
dangling portions. 

iii) FINITENESS t A solid must occupy finite portion 
of the space. 

iv) CLOSURE UNDER BOOLEAN OPERATION i Rigid motions or 
boolean operations when applied to solids must pro- 
duce other solids. Drilling, welding, assembly and 
design operations in general can thus be simulated 
in this way. 

v) FINITE DISCRIBABILITY t There must be some finite 
aspect of E^ (subset of 3-D Euclidean space, e.g. 
finite number of faces) to ensure that they are 
representable in computers. 
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1.3 Representation Schemes 

Traditional way of representation of the solid is by 
projection of solid, but this is an ambiguous representation 
unless one gives sufficient views of the same solid to 
remove ambiguities. The ambiguity arises because of the 
loss of information that occurs when a 3-D object is repre- 
sented on 2-D picture plane. 

Some of the unambiguous representations are* 

a) Pure Primitive Instancing Scheme ! Using this scheme 
one can represent only a limited range of the objects. 
Its drawback is that it cannot combine the created 
objects to model complicated objects and it is very 
difficult to compute properties of the solid repre- 
sented by this scheme. The advantage of this scheme 
is that it is concise and easy to use. It is like a 
language defined by its grammar. However, in this 
scheme, it is not possible to have explicit commands 
for boolean operations. 

b) Spatial Occupancy Enumeration ? Representation of 
solids in this scheme is essentially a list of cells 
occupied by a solid. The cells sometimes called 
’’Voxels*’ are cubes of fixed size which lie on a fixed 
spatial grid. Each cell may be represented by coor- 
dinates of a single point namely the centroid. It is 
mostly used to represent complex biological objects 
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encountered in computer aided tomography [9] • 

Octree encoding is a special case of spatial 
enumeration scheme. Octree is 3-D extension of 2-D 
quadtree encoding concept. In this scheme, a volume 
is divided into 8 octants till each part is either a 
homogeneous cell or contains a Voxel (Figure 1.1 ). 

c) Constructive Solid Geometry : Constructive solid 

geometry (C.S.G.) representation of the object is a 
binary tree whose nonterminal nodes represent operators 
which may be rigid motion like translation/rotation 

or regularised set operators like union intersection 
or difference. Every nonterminal node also represents 
a solid. The terminal nodes are the building blocks 
(primitive solids) such as cylinder, box, wedge etc. 

The basic idea is that complicated solids can be 
represented as various ordered additions and substrac- 
tions of some simpler solids or primitives. In this 
approach, the modelled object is stored as combination 
of data and logical procedures (boolean operation). 

This generally requires less storage but more compu- 
tation to reproduce the model. Figure 1.2 illustrates 
the C.S.G. representation [1]. 

d) Boundary Representation : In this scheme, a solid is 

represented by segmenting its boundary into finite 
number of bounded subsets called faces, edges, 






fig-1'2 constructive solid geometry 

REPRESENTATION 
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vertices. Boundary representation scheme is unambi- 
guous if the faces are represented unambiguously [2]. 

Following are the combinatorial conditions for 
true boundary representation: 

i) Each face must have atleast three edges, 

ii ) Each edge must have precisely two vertices, 

iii) Each edge must belong to even number of faces, 

iv) Each vertex in a face must belong precisely to 

the two of the face's edges. 

In boundary representation all the information is 
explicitly stored by defining the model geometry (faces, 
edges, vertices etc*). This requires more storage space but 
comparatively less computational effort to construct the 
image. This is helpful in modeling very complex geometries 
(Figure 1.3). 

e) Sweep Representation F 2] t In this scheme a solid 

object is represented by moving a set along a trajec- 
tory in space which sweeps a volume. There are two 
types of sweep operations. Translational sweep and 
Rotational sweep (Figure 1.4). Simple sweep represen- 
tations are easy to construct. But the domain of 
the objects they can describe contains only those 
objects that exhibit translational or rotational 
symmetry. Various generalised forms of sweeping are 
potentially capable of covering much large domains. 





9 



Translational sweep 





Generalised sweep producing undesired object 


FIG. 1*4 VARIOUS SWEEP REPRESENTATIONS 
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but they pose various mathematical and computational 
issues. Generalised sweeping sometimes creates 
unwanted objects with dangling edges or faces. 

1 .4 Statement of the Problem fs] 

One can create all the objects to be designed by 
manipulating a small library of fundamental solids (primi- 
tives) whose properties are well understood. 

Most of the solid modelers like GMSOLID, PADL, 

EUKLID, TIPS [l] use quadric solids such as cylinders, 
spheres, cones, alongwith boxes and wedges. 

These primitives cover a restricted range of solids 
and do not give the designer sufficient flexibility for 
modeling. 

There is, however, a new family of parametric 
objects called ' superquadrics' . These are generalizations of 
conic sections with the exponent 2 replaced by an appropriate 
positive number, n. These are flexible enough to represent 
large families of shapes which can be rounded, square etc. 
These can have different properties in different directions. 
The intent in the present woodc is to develop a geometric 
modeling scheme using * superquadrics ' . 

To understand the concept, consider the following 
equations. 
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X 

y 


a cos 



a 


0 p < 360 


( 1 . 1 ) 


Figure 1 .5 shows how N controls the shape of the curve. For 
N = 2, it is a circle and as N approaches oo the curve 
becomes scpjarish. 

This interesting behaviour was the source of motiva- 
tion for taking up the present work. In the present thesis, 
geometric modeling system based on constructive solid geometry 
scheme has been developed. The primitives used are super- 
quadric solids replacing the generally used quadric solids. 


1.5 Organisation of Present Work 

Chapter 2 gives brief summary of the Ray tracing 
theory, camera model for viewing, different primitives, 
coordinate systems, and the shading model. 

Chapter 3 deals with the boolean operations, 
algorithms for creating an object, boolean operation, Ray 
object intersection and data structure used for representing 
primitives and solids. 

Chapter 4 gives the basic structure of the program 
developed, different interactive features possible, proce- 
dures to be followed for creating/editing the modelled object, 
graphic displays and implementation details. 

Chapter 5 gives summary and scope for further work. 


r 
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CHAPTER 2 
RAY TRACING 


2.1 Ray Tracing Principle 

Ray tracing is essentially a photographic process 
in reverse. For each pixel (picture cell) element on the 
screen a light ray is cast into the scene. 

This ray acts as a sort of tracer. By tracing along 
the ray, ray-object intersections are calculated. These 
intersection points form the database which is used for mode- 
ling the objects. Modeling is achieved through boolean 
operations. Generation of shaded pictures or line drawings, 
calculation of various properties of the objects (such as 
volume, mass, moment of inertia) become a simple affair even 
for complicated objects. 

2.2 Camera Model f?] 

For viewing an object in any desired configuration, 
there are two equivalent ways of going about this* 

i) One can think of the view plane (in this case, the 

screen of the display device) to be fixed and object 
rotated for the desired configuration. 

ii ) By keeping the object fixed one can position the 
view plane according to the requirement (camera 
model). 
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For the present work second approach is used. This 
reduces the mathematical complexity involved in Ray-solid 
intersection procedure. 

Various viewing parameters used are: 

a) View Reference Point : This is the most important 

parameter in the camera model. It is a point in the 
scene which decides the portion of the scene that will 
be displayed on the screen. It acts as a sort of 
reference as all other parameters are taken with 
reference to this only. 

b) View Plane Normal Vector : View plane normal vector 

gives the direction perpendicular to the view plane 
directed towards the object. This means that the 
camera is always along the normal towards view refer- 
ence point. 

c) View Distance : It is the distance between the view 

reference point and the view plane measured along the 
view plane normal vector. 

d) View Up Direction : This parameter fixes the camera 

angle in the view plane. By changing this parameter 
one can imagine the camera pointing at the same object 
from the same direction but the picture will be 
tilted accordingly (Figure 2.2). 

By changing these parameters one can select different 
views of the object (Figure 2.1). 




View reference 
point 


View plane 


FIG. 2-1 DIFFERENT VIEWING PARAMETERS 
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2.3 Coordinate Systems and Transformations 

There are three different coordinate systems used 
for ray tracing technique. These arei 

i) Global Coordinate System : This is the most important 

coordinate system because in this coordinate system, 
user positions the camera to select a particular view, 
positions source of light for shading the object. 

User creates the object, performs boolean operations 
and modifies the object - all in this coordinate 
system only. 

ii ) Primitive’s Local Coordinate System : This coordinate 

system is primarily used to simplify the geometric 
calculations involved in ray-solid (primitive) inter- 
section. 

Each primitive has a local coordinate system. When- 
ever a user creates new primitive solid, its position and 
scaling in the global coo3:dinates are same as those in the 
local coordinates. But when a user moves or scales it then 
its position or size changes in the global coordinate system 
only. ' Local- to-global- transformation-matrix’ (Appendix A) 
is the link between the two coordinate systems. 

Ray-solid intersection can be found out in primitive's 
local coordinates and subsequently the intersection point can 
be transformed to global coordinates. 

iii ) Screen Coordinate System: Generally the screen or 
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the view plane is in the X-Y plane, the scene (object) 
is in the +Z half space and the viewer is in the -Z 
half space. 

The ray is generated in the screen coordinates. It 
is transformed to primitive' s local coordinates, through 
certain in-between transformations for the calculation of 
ray-solid intersection points. 

The block diagram for various transformations is 

given as 


Camera 
to global 


Global to object 




As parallel projection is used in this work, the 
direction of ray in the global coordinate is along view normal 
camera- to-global- transformation matrix is given in Appendix A. 
Global-to-obj ect- transformation matrix is given in Appendix A. 


2.4 Ray Primitive Intersection [5] 


A ray is a simple 3-D line in the camera model* It 
can be best defined in the parametric form as a point 
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P(XO, YO, ZO) and a direction vector V(DX, DY, DZ). In this 
form points on the line are ordered and accessed by a single 
parameter t. For every value of t the corresponding point 
on the line is given by 


X = XO + t DX 

y = YO + t DY (2.1 ) 

z = ZO + t DZ 

For transforming the ray from one coordinate system 
to another, only the point and vector transformations are 
essential. 

If [T] is 4 X 4 linear transformation matrix then 
P'(X', Y', Z', 1) *= P(X, Y, Z, 1) ♦ [t] (2.2) 


Similarly for the direction vector 


V'(DX', DY', DZ’, 1) = (DX, DY, DZ, 1) ♦ [t] 

(2.3) 


So, a line is simply transformed by 


(XO’, YO’, ZO’, 1)(DX’, DY’, DZ’, 1) = 

((XO, YO, ZO, 1) * [T])(DX, DY, DZ, l) ♦ [t] ) ) 

(2.4) 

From the above equations it is obvious that the line 
parameter ’t’ is same in both the coordinate systems. 
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This is very important because ray-primitive inter- 
section procedure finds only the parameter ' t' • For the 
calculation of intersection points solid or surface equations 
are not transformed, instead all the transformations are done 
on the ray. The effect of scaling, rotating and translating 
a solid is achieved by performing the above operations only 
on the ray. This simplifies the ray-primitive intersection 
calculation to great extent. 

The block diagram for transformation is given in 
Section 2.3. 

2.5 Primitive Definition [s] 

As this work is based on C.S.G, representation, there 
are certain basic building blocks or primitives used which 
satisfy all the criteria of the •’Abstract Solid" given in 
Section 1.2. 

The primitives used are: 

a) Box ; It is a unit cube in the positive octant with 
one corner at origin (Figure 2.3). 

The surface equations and the bound tests are 

X = 0 0:^Y, Z<1 

X = 1 0 < Y, Z ^ 1 

Y = 0 0:1X, Z;11 


Y 


1 


0 i: X, Z 1 
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Box 




FIG- 2-3 PRIMITIVES IN LOCAL COORDINATE 
CONFIGURATION 
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Z = 0 0<X, Y^1 

Z = 1 0 ^ X, Y 1 


To get the ray-box intersection we have to perform 
6 ray-plane tests. 

b) Wedge t The wedge is shown in the Figure 2.3. It has 

length, breadth and height equal to unity. 

The surface equations and the bound tests are 


X = 0 

X + Z = 1 

Y = 0 

Y = 1 

Z = 0 


0 ^ Y, Z < 1 
0 1 Y, Z ^ 1 

X > 0 
X + Z ^ 1 
Z 2. 0 

X 2 0 
X + Z < 1 
Z > 0 

0 1 X, Y :< 1 


To get ray-wedge intersection points we have to 
perform 5 ray-plane intersection tests. 

c) Cylinder ; It is a superquadric cylinder with its 
axis along the positive Z axis with one end at the 
origin (Figure 2.3). 

The surface equations with the bounds are 
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X^ + Y^ = 1 0 ^ ^ 

Z = 0 < 1 

Z = 1 1 

To get the intersection points two ray-plane inter- 
section tests and one ray- super quadric surface intersection 
test, need to be performed. 

d) Cone t It is a superquadric cone having its axis along 
the positive Z direction with the vertex at the origin 
(Figure 2.3). 

The surface equations alongwith the bounds are 

X^ + Y^-Z^ = 0 O^Zil 

Z = 1 X^ + Y^ £ 1 

For the ray-cone intersection points we have to 
perform one ray-plane and one ray— sup er quadric surface inter- 
section test. 

e) Sphere t It is a superquadric sphere having its radius 

unity and the centre at the origin (Figure 2*3). 

The surface equation is 

+ = 1 

To get ray-sphere intersection we have to perform one 
ray— sup er quadric surface intersection test. 
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2.6 Rendering 

Rendering is essentially producing the realistic 
pictures or images. 

In a simple shading model (Lambertian reflection 
model) [s] a source of light is positioned with respect to the 
object. For every visible surface point on the object, light 
intensity is calculated by calculating the angle between the 
normal and the line joining the surface point and the light 
source. 

Objects rendered with a simple Lambertian diffuse 
reflection illumination model appear to have dull matte 
surface because a point light source is assumed. Objects 
that receive no light directly from the source appear black. 
However in real scene, objects also receive light scattered 
back to th^ from surrounding. 

An empirical shading model given by Bui-Tong- Phong [s] 
is implemented in this work. This model is explained as 
follows (Figure 2.4). 

If I is the intensity of the surface as viewed by the 
observer then 

I = ♦ (L . N) + Kg (V . nf 

where 

t diffused reflection constant 
N : unit surface normal 
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K : specularity constant 

o 

A 

L : unit light vector 

V : L - S (S - unit sight vector) 

E : constant which determines the spatial distri- 
bution characteristics of metallic or non- 
metallic surfaces. 
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CHAPTER 3 

BOOLEAN OPERATIONS AND ALGORITHMS 

3.1 In Out Classification Fs] 

For a given prlniitive, ray-primitive intersection has 
four possible outcomes t 

1 ) Ray missing the primitive 

2) Ray touching the primitive at one point 

3) Ray enters and exits through primitive at two different 
points 

4) Ray lies on face of a primitive. 

For cases (l ) and (2), a ray is classified as comp- 
letely out and so no intersection points are noted. In cases 
(3) and (4), two intersection points are noted and the ray is 
divided into 3 segments: out-in-out. This classification 
is very important for carrying out the boolean operations. 

3.2 Data Structure 

As mentioned earlier the data structure representing 
the constructive solid geometry scheme is an inverted binary 
tree. 

The terminal nodes of the tree are the different 
primitives used. The internal nodes are the composite solids 
formed by different boolean operations. 

The data structure used to represent a primitive 
solid has a parameter attribute. This stores the information 
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such as superquadric degree, rotation/translation of the 
primitive about the axes, different dimensions such as length, 
breadth, height, major and minor radii depending on primitive 
type. 

It also has world-to-obj ect and object-to-world 
transformation matrices as attributes. These are used for 
representing the primitive in global coordinates or transfor- 
ming the ray in primitive' s local coordinate system. 

The composite solid has operation attribute which 
specifies the boolean operation UNION (+), DIFFERENCE (-) or 
INTERSECTION (<^). It has left and right pointers for pointing 
to left subtree and right subtree. The last attribute is leaf 
number. If this is zero then the node is a terminal node. 

All the objects created Including the primitive 
solids are accessed by pointer which points to a particular 
node of the tree depending upon the name attribute. 

3.3 Boolean Operations 

Solid primitive and boolean operation comprise a 
natural and powerful design mode. For example, when a designer 
wants to drill a hole in a part what could be more natural 
than to place a cylinder in the desired position of the hole 

and then substract it from the part. 

In a C.S.G. representation solid objects are modelled 
as composition of primitive solids using boolean operations 
namely union, intersection and difference. 
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The ray solid intersection procedure starts at the 
top of the tree, recursively descends to the bottom classi- 
fying the ray with respect to primitive solid and then 
returns up the tree combining the classification of the left 
and right subtrees (Figure 3.1). In Figure 3.1 dotted portion 


of the ray is 

inside the solid 

and solid 

portion is outside 

the solid. 




The boolean operations 

are simple 

set operations. 

The rules are 

given below: 



Operator 

Left 

Right 

Composite 

Union 

IN 

IN 

IN 


IN 

OUT 

IN 


OUT 

IN 

IN 


OUT 

OUT 

OUT 

Intersection 

IN 

IN 

IN 


IN 

OUT 

OUT 


OUT 

IN 

OUT 


OUT 

OUT 

OUT 

Difference 

IN 

IN 

OUT 


IN 

OUT 

IN 


OUT 

IN 

OUT 


OUT 

OUT 

OUT 


To explain these operations more clearly consider the 
example as shown in Figure 3.2. Solid line represent ray 
portion inside the solid and dotted indicates that ray is out- 
side the solid. 
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3.4 Algorithm for Creating the Object 

This procedure as the name suggests when invoked 
creates a new object. Different cases for creating the new 
objects are: 

1 ) Two primitives can be combined to form a new object 

2) One object (already created) and one primitive can be 
combined to form a new object 

3) Two objects (already created) can be combined to form 
a new object 

4) Already existing object can be modified. 

These features give the user the flexibility needed 
to create new objects. 

The algorithm is presented below; 

a] Read the object name from TTY and store the object in 
a stack; get the corresponding stack value. 

b] Assign the stack value to a variable LHS. 

c] If LHS is less than zero then 

i ) prompt error message 

Attempt to redefine primitive 
ii ) exit from procedure. 

d] Read a character from TTY 

e] If character is not equal to '=’ then 

i ) prompt error message 
Equal to *=' expected 


ii) exit from procedure 



Read the object name from TTY and pop the object stack 
to get the corresponding stack value. 

Assign the stack value to a variable RH1 . 

If RH1 is equal to zero then 

i) prompt message 
Undefined objecti try again 

ii) exit from procedure, 
otherwise 

if RHI is less than zero then 

i) assign the parameters if any to RH1 

ii) create a node with pointer to left 

iii) assign the node the corresponding leaf number 
otherwise 

i) make a copy of the objects of the left subtree 
ii) if parameters then modify the object. 

Read a character from TTY. 

If character is equal to then 
If RH1 is less than zero then 
i ) create a new node 

ii ) assign the corresponding leaf number to 
this node 

iii ) store the obj ect in the obj ect stack 
otherwise 

i) store the modified object in the object list 
ii) set the pointer to indicate .to left subtree. 



35 


k] Read the boolean operation from TTY. 

l] Read the object name and get the corresponding stack 
value for this object. 

m] Assign stack value to variable RH2. 

n] If RH2 is equal to zero then 

i) pTOmpt message 
Undefined object! try again 

ii ) exit from procedure 
otherwise 

If RH2 is less than zero then 

i) assign the parameters if any to RH2 

ii ) create a node with pointer pointing to right 

iii) assign the node corresponding leaf nuiaber 
otherwise 

i) make a copy of the objects in the right subtree 

ii) if parameters then modify the object. 

o] Create a node, assign the left pointer to the left sub- 
tree, right pointer to the right subtree. 

p] Make an entry of the node in the object stack. 

3.5 Algorithm for Boolean Operations 

This procedure combines the left subtree intersection 
points and right subtree intersection points according to the 
boolean operation. 
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The parameters required are boolean operation, inter- 
section points in left subtree and right subtree* The output 
is the combined intersection points depending on the boolean 
operation. 

Let LST1 and LST2 be the two stacks containing left 
subtree and right subtree intersection points respectively 
and LST is the resultant stack after boolean operation. 

The algorithm is then as follows: 

a ] If left stack is empty then 
Case boolean operation of 

i) Union : LST *■ right stack (LST2) 

ii) Intersection : LST null stack 

iii) Difference : LST null stack 

otherwise 

b] If right stack is empty then 
Case boolean operation of 

i) Union : LST •- left stack (LST1 ) 

ii ) Intersection : LST •• null stack 

iii) Difference t LST .* null stack 

otherwise 

c] Pop the element out of right and left stack 

d] While left stack or right stack is not empty do 

I] If left and right intersection points are same then 
If left subtree stack is empty then 
i) assign right stack value to temporary list 
ii ) pop the right stack for next value 
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otherwise 

If right subtree stack is empty then 
i) assign left stack value to temporary list 
ii ) pop the left stack for next value 
otherwise 

i) assign left stack value to temporary list 
ii) pop right and left stack for next value 
otherwise 

II] If left stack value is less than right stack value 
then 

i) assign left stack value to temporary list 
ii ) pop left stack for next value 
else 

i) assign right stack value to temporary list 
ii ) pop right stack for next value 
III] Store the temporary list value in the resultant 
stack (LST) depending on boolean operation 

e] If left stack or right stack is not empty then go to 
step (d). 


3.6 Ray-Object Intersection Algorithm 

This procedure starts at the top of the tree, descends 
downwards to the terminal nodes, calculates the intersection points 


and returns up the tree combining left and right subtrees 
The algorithm is as follows: 
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a] For each pixel in the window do 

b] Set the view direction 

cj Initialise the intersection point stack 
d] If pointer is not pointing to a primitive then 

i ) get left subtree intersection stack 

ii ) get right subtree intersection stack 

iii ) merge the two stacks according to boolean operation, 
else 

i) transform the ray to primitive's local coordinate 
syst®n by primitive's world- to-obj ect transfor- 
mation 

ii ) do case primitive type of 
BOX : Solve 6 ray-plane intersection test 

Wedge : solve 5 ray-plane intersection test 

CYLINDER : Solve 2 ray-plane intersection and one ray 

superquadric surface intersection test 
CONE : Solve one ray-plane and one ray-superqua- 

dric surface intersection test 
SPHERE : Solve one ray- superquadric intersection test. 

3.7 Bairstow* s Method [b] 

For ray- super quadric surface intersection we have to 
solve a polynomial for the intersection points. For solving 
the polynomial Bairstow' s method is iraplesented. This proce- 
dure breaks the polynomial in quadratic parts and then calcu- 
lates the roots of the quadratic part if any. 
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The algorithm is as follows. 

Let a^ contains the coefficients of polynomial and 
n is the total number of terms in the polynomial (1+ degree 
of the polynomial). 

‘ W2 

b] Set 

b^ * 1.0; b^ * 32 - U 
Ci ^1.0; C 2 ^ b2 - U 

c] Calculate bj^ aj^ - aj^ ^ ♦ U - aj ^_2 * V 

d] Calculate Cj^ \ ~ * U - Cj ^_2 '**’ ^ 

e] Calculate Deno *■ c^_^ * ®n-3 “ ^n-2 

- '=nV2>''°®“'’ 

ERROR - ABS(DU + DV) 

f] If Error is greater than or equal to converging limit 
then go to (b). 

g] Calculate root of quadratic part 

h] If the remaining part of the polynomial is not quadratic 
or linear then go to (a). 

i] Calculate the root of the quadratic part or linear part 
as the case may be. 


Exit 
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CHAPTER 4 
SOFTWARE development 

4. 1 Implementation Details 

The present software has been developed in PASCAL 
language and Implemented on ND 560/CX system* For graphic 
displays Tektronic 4109 and 4107 raster graphic terminals are 
used. For this software various local commands of the term- 
inal are interfaced for different graphic procedures. 

4.2 Program Structure 

Figure 4.1 gives block diagram of the software deve- 
loped. Initiate data procedure initiates certain data such as 
it sets number of objects created to zero, selects a default 
window etc. 

Initiate device procedure activates the terminal, 
clears the graphic displays present on the screen. 

The following section give details of the various 
interactive features. 

4. 3 Interactive Features 

The basic requirement of any geometric modeling system 
is that it should be highly interactive and user's friendly. 
The user should have maximum possible freedom for modeling as 
well as viewing the object. Keeping this as a guide line. 
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certain commands are developed which are simple and highly 
interactive. 

When executed the program displays the message and 
the prompt 

TYPE HELP (FOR HELP) 

'=>' This prompt is very important because it indicates that 
the program is waiting for the next command. 'Whenever a . 
particular command is over, '=^' is displayed. 

4.3.1 Create ; This command is most important. With this 
command only, the user creates different objects by combining 
two already existing objects or by combining a object and a 
primitive or by combining two primitives. User can modify 
an existing object by translating/rotating it. 

The primitives along with the associated parameter are: 
i) BOX(L, B, H, T, P, W, X, Y, Z, D, S, E, C): 

L = Dimension along X direction (an integer) 

Default = 100 

B = Dimension along Y direction (an integer) 

Default = 100 

H = Dimension along Z direction (an integer) 

Default = lOO 

T = Rotation about Z axis in degrees (an integer) 
Default = 0 

P = Rotation about Y axis in degrees (an integer) 


Default = 0 
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W = Rotation about X axis in degrees (an integer) 

Default = 0 

X = Translation along X axis (an integer) 

Default t= 0 

Y = Translation along Y axis (an integer) 

Default = 0 

Z = Translation along Z axis (an integer) 

Default = 0 

D = Object light diffusivity (an integer) 

Default = 70 

S = Object light specularity (an integer) 

Default = 30 

E = Object specular exponent (an integer) 

Default = 2 

C = Colour for the line for drawing the wire frame of 
obj ect 

Default = 15 

ii) WEDGE(L, B, H, T, P, W, X, Y, Z, D, S, E, C) 

All the parameters are same as discussed for BOX. 

iii) CYLINDER(N, A, B, H, T, P, W, X, Y, Z, D, S, E, C) 

N = Superquadric degree (an integer) 

Default = 2 

A = Major axis of the cross section (an integer) 

B = Minor axis of the cross section (an integer) 

H = Height of the cylinder 

Rest of the parameters are same as given for the BOX, 
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iv) C0NE(N, a, B, H, T, P, W, X, Y, Z, D, S, E, C) 

All the parameters are same as discussed for 
cylinder, 

v) SPHERE(N, R, T, P, W, X, Y, Z, D, S, E, C) 

R = Radius of sphere (an integer) 

Default =100 

Rest of the parameters are same as given for, 
cylinder. 

The characters used for boolean operations are 

'+' for Union 
*-* for Difference 

for Intersection. 

Different cases elaborating how to give the input for 
create routine are given below: 

a) Input for combining two premitives 

=> CREATE PIPE = CYLINDER (A = 40, B = 40, H = 100) - 

CYLINDER (A =30, B = 30, H = 100) 

By this We are creating an object with name attribute 
as pipe. This consists of two cylinders (parameters are 
listed) with the second substracted from the first. Those 
parameters whose values are not read are taken as default. 

b) Input for modifying an existing object 

A typical command is 
=> CREATE A. = B(T, P, W, X, Y, Z); 

Here B is an already existing object which we want to 
modify, A stores the modified definition of B, The 



45 


parameters which can be modified are listed. 

c) Input for combining two already existing objects 

=> CREATE A = B(T, P, W, X, Y, Z) + C(T, P, W, X, Y, Z) 

Here B and C are already existing objects. A stores 
the combination of B and C. Any of the parameters listed 
can be changed according to requirement. 

4.3.2 L^: This command gives all the relevant information 

regarding an object. It gives the list of all the primitives 
with their parameters and associated boolean operations which 
are used to model the solid. 

A typical command is 

=> LIST OBJECT-N/ME 

where OBJECT-NAME in general is the name of the object created. 
For illustration consider the example given in Section 4.3.1. 

For listing the command is 
=> LIST PIPE 

The response is as follows: 

CYLINDER (2)(0, 0, 0)(0, 0, 0)(40.00, 40.00, 100.00)(2, 70, 30) 

CYLINDER (2)(0, 0, 0)(0, 0, 0)(30.00, 30.00, 100.00)(2, 70, 30) 
This feedback is very essential while modeling a complex 

obj ect. 

4.3.3 Draw ; This command draws the wireframe picture of the 
object. It gives the user graphical feedback regarding the 
problem during modeling. User can select different colours for 
different primitives for clear perception. 
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The command is 

=i> DRAW OBJECT 

This will draw all the primitives associated with 

OBJECT. 


4.3.4 Di_stancet This command specifies the position of 
view plane with respect to an object or scene. The distance 
should always be a positive integer as according to the conven- 
tion followed, the scene is always in front of viewer. A typi- 
cal command is 

=> DISTANCE 20 

This sets the view plane along the view plane normal 
at a distance of 20 units. 

4.3.5 Vi ewup t This command specifies the viewup vector 

reguired for camera setting. The eluents of the vector are 
real numbers. A typical command is 

^ VIEWUP 0 0 1 

This sets the camera angle in the view plane. 

4.3.6 Normal ! This command specifies the view normal direc- 
tion. The vector elements to be keyed in must be integers. A 
typical command is 

=> NORMAL -1 -1 -1 

This sets the viewer’;s direction in the global coordi- 
nate system. 

4 . 3.7 Position : This command sets the view reference point. 

All the coordinates to be specified are integers. 



47 


A typical command is 

=i> POSITION X coordinate Y coordinate Z coordinate 
The default values are 000. 

4.3.8 Axes t This command draws the global coordinate system 
when invoked. This command should be invoked only after 
setting the desired camera parameters. The command is 

^ AXES 

4.3.9 Window: This command reads the window parameters which 

are in screen coordinate system. This helps to select different 
windows according to this requirements. A typical command is 

=> WINDOW X min Y min X max Y max 

All the window parameters are positive integers. 

4.3.10 Light : This command sets the light source in the 

world coordinate system. The typical command is 

LIGHT X coordinate Y coordinate Z coordinate 

4.3.11 Remove : This command removes the light source existing 

thus giving the user the flexibility to again place the source 
of light at desired location. The command is 

^ REMOVE 1 

4.3.12 Colour : This command reads the hue value and accor- 

dingly sets the colour indices. These, indices are subsequently 
used for shading the object. Typical command is 


=> COLOUR VALUE 
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B^.sterize: This is one of the important commands. 

When invoked this produces the solid model of the object in 
the window selected. All the operations like setting the 
camera, creating the desired object, selecting the proper 
window, setting the light source should be done before this 
command. A typical example is 

=*■ RASTERIZE OBJECT 

If the object is not defined or the name attribute 
keyed in is incorrect then a error is prompted 
UNDEFINED OBJECT! TRY AGAIN. 

4.3.14 Help t As the name suggests, when this is invoked, a 
list of all the interactive commands available is provided to 
the user. To invoke help 

HELP 

4.3.15 Exit : To stop the program this command is to be 

invoked. It allows the control to come out of the interactive 
repeat loop and stops the entire program. The command is 

=> EXIT 

4.4 Sample Displays 

Figure 4.2 shows the picture of Box and Wedge. Figure 
4.3 shows the picture of superquadric cylinder with superqua- 
dric degree 4 and general quadric cylinder. Figure 4.4 shows 
the picture of superquadric cone with superquadric degree 4 
and general cone. Figure 4.5 shows the picture of superquadric 




Pigure 4.2 



Figure 4.3 
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Figure 4*5 
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Figure 4.6 



Figure 4*7 
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Figure 4.8 
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sphere with superquadric degree 4 and general sphere. All 

above are the basic building blocks (primitives) incorporated 
in the present work, 

4.4.1 U nion t Figure 4.6 illustrates the union operation. 

The basic primitives used are box and superquadric cylinder. 

With reference to Section 4.3, for producing this object the 
input was given in the following manners 

=> DISTANCE 20 
VIEWUP 0 0 1 
NORMAL -1 -1 -1 
LIGHT 100 -200 -400 

^ CREATE OBJECT1 = B0X(L = 60, B = 60, H = 10); 

=t» CREATE OBJECT2 = CYLINDER(N = 4, A = 20, B = 20, H = 60, 

Z = —30, X = 30, Y = 30); 

=i> CREATE 

si> WINDOW 100 100 200 200 

=> RASTER OBJECT1 
=> WINDOW 300 100 400 200 

=?► RASTER OBJECT2 
=> WINDOW 200 200 300 300 

=> RASTER OBJECTS 

4.4.2 Difference s Figure 4.7 illustrates the difference oper- 
ation. The basic primitives used are box and superquadric 
cylinder. 

With reference to Section 4.3, for producing this 
picture input was given in following manners 


54 


=> DISTANCE 20 
=> VIEWUP 0 0 1 
NORMAL -1 -1-1 
=4> LIGHT 100 -200 -400 

=> CREATE OBJECT1 = BOX(L = 60, B = 60, H = 10); 

=> CREATE OBJECT2 = CYLINDER(N = 4, A x= 20, B = 20, H = 60, 

Z = -30, X = 30, Y = 30); 

=> CREATE OBJECT3 = OBJECT 1 ( ) - OBJECT2( ) 

=> WINDOW 100 100 200 200 

=* RASTER OBJECT1 

WINDOW 300 100 400 200 

=> RASTER 0BJECT2 

=> WINDOW 200 200 300 300 

RASTER OBJECTS 

4.4.3 Intersection t Figure 4.8 illustrates the intersection 
operation. The basic primitives are the same i.e. box and 
superguadric cylinder. With reference to Section 4.3 for 
producing this object the input was given in this manner; 

=> DISTANCE 20 
VIEWUP 0 0 1 
=> NORMAL -1 -1 -1 
=> LIGHT 100 -200 -400 

=> CREATE OBJECT1 = B0X(L =60, B = 60, H = 10); 

=S> CREATE 0BJECT2 = CYLINDER(N = 4, A = 20, B = 20, H = 60, 

Z = -30, X = 30, Y = 30); 

=> CREATE OBJECTS = OBJECTI ( ) * 0BJECT2( ) 

=> WINDOW 100 100 200 200 


=> RASTER OBJECT 1 

WINDOW 300 100 400 200 
RASTER OBJECT2 

=> WINDOW 200 200 300 300 


RASTER OBJECTS 
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CHAPTER 5 
SUMMARY 


5.1 Conclusions 

A geometric modeler based on constructive solid geo- 
metry representation scheme using superquadric solids as primi- 
tives has been implemented and analysed. 

Many times in engineering practice it is observed that 
the workpieces do not have sharp edges* These classes of work- 
pieces are difficult to model using generally used prism and 
quadric solid primitives. It is observed that the use of the 
superquadric solids gives the user greater flexibility in 
modeling such objects. These superquadric solids cover a wide 
range of objects which are very useful in modeling different 
complex objects. It is possible to modify the basic primitive 
itself by changing the superquadric degree according to the 
requirement. 

5.2 Suggestions for Further Work 

a) Ray tracing in general qualifies as a brute force 

method for solving problems. Memory and CPU usage is directly 
proportional to the scene complexity. 

Most of the efforts are concentrated on calculating 
ray-solid intersection points. Even though there is no object 
in the path of the ray still it visits each node of the tree 
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and performs all the intersection tests. By using a minimum 
bounding box around a solid in the composition tree, the 
exhaustive search for a ray-solid intersection then resembles 
an efficient binary search. This improves the efficiency of 
the ray-tracing procedure to great extent. So bounding box 
technique can be incorporated with present work to improve the 
efficiency of ray-solid intersection procedure. 

b) Various application subsystems for calculating mass, 
moment of inertia volume which are essential for a solid 
modeler can be interfaced with the present work. This would 
result in a complete solid modeler. 

c) In the present approach, the complex object is created 
by combining different primitives through certain operations. 
The same approach can be applied in reverse manner also. 

Given a complex object composition one can get back the basic 
primitives used to model the object. This will be very helpful 
for deassembling operation of the mechanical components. 
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APPENDIX A 

TRANSFORMATION MATRICES 

A1 Object-to-Global Transformation Matrix 


Let SX, SY and SZ be- the scaling factor in X, Y and 
Z direction. 0, 0 and w are the angle of rotations in radians 
about Z, Y, X axis. X, Y, Z be the translations along X, Y, Z 
axis then the transformation matrix [T] is given by 


cos0.cos0.SX SX(sin0.cosw.— 

cos0.sin0.sinw) 

-sin0.cos0.SY SY(cos0.cosw + 

sin0.sin0.sinw) 


-SZ.sin0 


-SZ.cos0.sinw 


SX(sln0.sinw + 

CO s0 . sin0 . CO sw ) 

SY(cos0.sinw - 
sin0.sin0.cosw) 

SZ.cos0.cosw 


0 

1 


The inverse of the matrix [T] gives global- to-object 
transformation matrix. 


A2 Camera- to-World Transfoirmation Matrix 

The purpose of this transformation is to establish 
relationship between camera and world coordinate system. 

Let DXN, DYN, DZN be the component of view normal 
vector. DIS be the viewer' s distance. XUP, YUP, ZUP be the 
component of view up vector. XR, YR, ZR be the coordinates 
of view reference point. Then the transformation matrix 
is given by 
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M 


where 


[T] = 


[R^] = 


where 


1 

0 

0 




0 

1 

0 


0 

0 

1 


-(XR+DXN^^DIS) - (YR+DYN*DIS ) -(ZR+DZN*DIS ) 


1 

0 

0 

0 


0 

-DZN/V 

DYN/V 

0 


0 

-DYN/V 

-DZN/V 

0 


V = (DYN^ + DZN^)^'^^ 


0 

0 

0 

1 


0 

0 

0 

1 


[Ry] = 


[\] 


V 

0 

DXN 

0 

YUP/RUP 

-XUP/RUP 

0 

0 

RUP : 


0 

1 

0 

0 

XUP/RUP 

YUP/RUP 

0 

0 


-DXN 

0 

V 

0 

0 

0 

1 

0 


(XUP^ + YUP^)"*/^ 


0 

0 

0 

1 

0 

0 

0 

1 



